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Searches for lepton-flavor-violating decays of a r lepton to a lighter mass lepton and a photon 
have been performed with the entire dataset of (963 ± 7) x 10 6 r decays collected by the BABAR 
detector near the T(4S), Y(3S) and T(2S) resonances. The searches yield no evidence of signals 
and we set upper limits on the branching fractions of B(t ± — > 6*7) < 3.3 x 10~ s and B(t ± — > n ±n i) 
< 4.4 x 10~ s at 90% confidence level. 
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Amongst all the possible lepton-flavor- violating r pro- 
cesses, r ± — > £ ± 7 (where 1 = e,fx) is predicted to be the 
dominant decay mode in a wide variety of new physics 
scenarios, with rates close to current experimental lim- 
its. Despite the existence of neutrino oscillations [H, such 
decays are predicted to have unobservably low rates 
in the Standard Model (SM). Thus, an observation of 
charged lepton flavor violation would be an unambiguous 
signature of new physics, while improvements on existing 
limits will constrain many models. As the relationships 
between — > e ±r y, — > e ± j and — > /i ± 7 decays 
are model-dependent, searches for both r modes provide 
independent information, even in the light of the small 
limit of B(/i + e+7) < 1.2 x lCT 11 at 90% confidence 
level (C.L.) @|. 

Presently, the most stringent limits are B{t ± — > e 7) 
< 1.1 x 10" 7 Q and B{t± -» ^7) < 4.5 x 10~ 8 @ at 
90% C.L., using 232.2 fb _1 and 535 ftT 1 of e+e~ anni- 
hilation data collected near the T(45 l ) resonance by the 
BABAR and Belle experiments, respectively. This paper 
reports the final result on these processes from BaBar. It 
utilizes the entire dataset recorded by the BABAR detector 
at the SLAC PEP-II e + e~ storage rings, corresponding 
to a luminosity of 425.5 fb"\ 28.0 ftT 1 and 13.6 ftT 1 
recorded at the T(iS), T(3S) and T(2S) resonances, and 
44.4 fb"\ 2.6 ftT 1 and 1.4 ftT 1 recorded at 40MeV, 
30 MeV and 30 MeV below the resonances, respectively. 

For the bulk of the data sample at the T(45) res- 
onance, the cross-section cr e + e -_j. r + T - = (0.919±0.003) 
nb [f|, determined to high precision using the KK Monte 
Carlo (MC) generator [7j, receives negligible contribu- 
tion from T(4:S) due to its large decay width. But, for 
the remaining data at the T(3S) and T(2S) resonances, 
the r-pair cross-section receives additional contributions 
of B(T —> t + t~) ss 2%, which are known only at the 
10% level [U . Including a systematic uncertainty of 0.6% 
on the luminosity determination, this gives a total of 
N T = (963 ± 7) x 10 6 r decays. 

The BABAR detector is described elsewhere 
Charged particles are reconstructed as tracks with a 5 
layer silicon vertex tracker and a 40 layer drift cham- 
ber inside a 1.5 T solenoidal magnet. A CsI(Tl) elec- 
tromagnetic calorimeter is used to identify electrons and 
photons. A ring-imaging Cherenkov detector is used to 
identify charged pions and kaons. The flux return of the 
solenoid, instrumented with resistive plate chambers and 
limited streamer tubes, is used to identify muons. 

The signal is characterized by a ^7 pair with an in- 
variant mass and total energy in the center-of-mass (CM) 
frame (E™ 1 ) close to m T = 1.777 GeV/c 2 8] and y/s/2, 
respectively. The event must also contain another t de- 
cay, reconstructed as decaying to one or three tracks. 

The dominant irreducible background comes from r- 
pair events containing hard photon radiation and one 
of the t leptons decaying to a charged lepton. The re- 

e ± 7 and ^ — » /^ ± 7 de- 



maining backgrounds for r 



cays arise from the relevant radiative processes, e + e~ — > 
e + e~7 and e + e~ — > /i + /i~7, and from hadronic r decays 
where a pion is misidentified as the electron or muon. 

Signal events are simulated using KK and TAU0LA [l(| 
with measured r branching fractions 0]. The 
and t + t~ background processes are generated using KK 
and TAU0LA, while the qq processes are generated using 
JETSET [Tl| and EVTGEN[12|]. Radiative corrections for all 
processes are simulated using PHOTOS [ljj. The Bhabha 
background is studied using events with two identified 
electrons in the data. The two-photon background has 
been studied and found to be negligible. The detector 
response to generated particles is simulated using the 
GEANT4 package 14j. MC events are used to optimize 
the selection criteria and estimate the systematic uncer- 
tainties on the efficiency, while the background rates are 
estimated directly from data. 

Events with two or four well reconstructed tracks and 
zero total charge are selected, where no track pair is 
consistent with being a photon conversion in the de- 
tector material. Each event is divided into hemispheres 
( "signal-" and "tag-" sides) in the CM frame by a plane 
perpendicular to the thrust axis calculated using all re- 
constructed charged and neutral particles [15j. 

The signal-side hemisphere must contain one photon 
with CM energy E^ M greater than 1 GeV, and no other 
photon with energy greater than 100 MeV in the labora- 
tory frame. The signal side must contain one track within 
the calorimeter acceptance with momentum in the CM 
frame less than 0.77^/2. This track must be identified 
as an electron or a muon for the t — > e ± 7 or r — » /z 7 
search. The electron selectors have an efficiency of 96% 
within the fiducial coverage. For reliable muon identi- 
fication, the track momentum is required to be greater 
than 0.7 GeV/c in the laboratory frame, above which the 
selection efficiency is 83%. 

In the rest-frame of the , the ^ ± and the 7 are pro- 
duced back-to-back. When boosted to the CM frame, 
kinematic considerations of two-body decays require 
there to be a minimum opening angle between them. The 
cosine of the opening angle, cos^ 7 , between signal-track 
and signal-photon is required to be less than 0.786. 

The tag-side hemisphere is expected to contain a SM 
r decay. A tag-side hemisphere containing a single track 
is classified as e-tag, /i-tag, or 7r-tag if the total photon 
CM energy in the hemisphere is less than 200 MeV and 
the track is exclusively identified as an electron (e-tag), 
as a muon (/i-tag), or as neither (7r-tag). Events with 
the tag-side track failing both the lepton selectors are 
classified as p-tag if they contain at least one ir° candidate 
reconstructed from a pair of photons with invariant mass 
between 90 and 165 MeV/c 2 . If the tag-side hemisphere 
contains three charged tracks, all of which fail the lepton 
identification, it is classified as a 3/i-tag. 

The definitions of the tag-side modes are designed to 
minimize the residual backgrounds from radiative QED 
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processes. For the r — > e ± 7 search, very loose electron 
selection criteria are applied for the e-tag sample. Thus, 
the remaining tags which fail these very loose electron cri- 
teria have small Bhabha contamination. The e-tag events 
are used as the control sample to model the Bhabha back- 
ground characteristics, and are removed from the final 
sample of events in the r — > e ± 7 search. Similarly, 
for the — > /i ± 7 search, very loose muon criteria are 
applied for the p-t&g, on which stricter kinematic require- 
ments are later imposed with tolerable loss in signal effi- 
ciency. The other tags are required to fail these very loose 
muon criteria, thereby reducing di-muon backgrounds. 

To suppress non-r backgrounds with missing momen- 
tum along the beam direction due to initial and final 
state photon radiation, we require that the polar angle 
m i ss of the missing momentum be inside the detector 
acceptance, i.e. —0.76 < cos# m i SS < 0.92. 

The total CM momentum of all tracks and photon 
candidates on the tag-side is required to be less than 
0.77y/s/2 for e-, p,-, 7r-tags and less than 0.9-^/2 for p- 
and 3/i-tags. The tag-side pseudomass is required to 
be less than 0.5GeV/c 2 for e-, p-, 7r-tags and less than 
1.777 GeV/c 2 for p- and 3/i-tags. 

The mass squared m 2 of the missing particles on the 
tag side is calculated using the tag-side tracks and pho- 
ton candidates and assuming that in the CM frame, 
the tag-side r momentum is opposite that of the sig- 
nal r and that its energy is \/s/2. To reduce back- 
grounds, we require to 2 > — 0.25GeV 2 /c 4 for e- and 
/i-tags, |m 2 | < 0.25GeV 2 /c 4 for 7T- and 3/i-tags, and 
|m 2 | < 0.50GeV 2 /c 4 for p-tags. 

For radiative Bhabha and di-muon events, the ex- 
pected photon energy in the CM frame (E^ M ) cxp is 

■ a r^arr^-nr i a m > where n— 9i and n— 6> 2 are the an- 
gles the photon momentum makes with the signal-track 
and the total observed tag-side momentum, respectively. 
Also, for such events, we expect the cosine of the open- 
ing angle, # r ecoii, between the signal-track and the to- 
tal observed tag-side momentum in the reference frame 
obtained by removing the signal photon from the CM 
frame to peak at -1. To suppress these backgrounds, we 
remove events having reconstructed photon energy con- 
sistent with the expected value, i.e. \E^ M - (£^ M ) cxp | = 
|A£J 7 | < 0.125^ and cos0 rcC oii < —0.975 in e- and p- 
tags for the r — > // 7 search. No such criteria are 
necessary for the r — > e ± j search according to the op- 
timization procedure. 

To further suppress the remaining backgrounds, neural 
net (NN) based discriminators are employed separately 
for each tag and for each dataset taken at values of ^/s 
near the T(45), T(3S) and T(2S) resonances. Six ob- 
servables are used as input to the NN: the total tag-side 



7/ V s ) COS "recoil, 



momentum divided by \fs/2, m 2 ,, AE n 
cos 0£j , and the transverse component of missing momen- 
tum relative to the collision axis. The NN based discrim- 



inators improve the signal to background ratios for the 
two searches by factors of 1.4 and 1.3, respectively. 

Signal decays are identified by two kinematic variables: 
the energy difference AE — E™ — \J~sj2 and the beam- 
energy constrained r mass (raEc), obtained from a kine- 
matic fit after requiring the CM r energy to be y/s/2 and 
after assigning the origin of the 7 candidate to the point 
of closest approach of the signal lepton track to the e + e~ 
collision axis. The distributions of these two variables 
have a small correlation arising from initial- and final- 
state radiation. For signal MC events, the 77t.ec and AE 
distributions are centered at m T and small negative val- 
ues, respectively, where the shifts from zero for the latter 
are due to radiation and photon energy reconstruction 
effects. The mean and standard deviations of the tbec 
and AE distributions for the reconstructed signal MC 
events are presented in Table Q] The data events falling 
within a 3cr ellipse in the ttzec vs. AE plane, centered 
around the reconstructed peak positions as obtained us- 
ing signal MC, are not examined until all optimization 
and systematic studies have been completed. The selec- 
tions are optimized to yield the smallest expected upper 
limits 17J for observing events inside a 2a signal ellipse 
under background-only hypotheses. 

The distributions of events in ttiec vs. AE are shown 
in Fig. [1] To study signal-like events, a Grand Signal 
Box (GSB) is defined as m EC € [1.55,2.05] GeV/c 2 and 
AE e [-1.0,0.5] GcV. Outside the blinded 3a ellipse, 
1389 data events survive in the GSB for the r — > e ± 7 
channel, and 2053 data events survive for the t —> // 7 
channel. These agree to within 2.4% and 1.7% with the 
numbers of background MC events observed. The signal- 
track arises from a real electron or muon in 96% and 82% 
of the background MC events for the two searches. 

A Fit Box (FB) region is defined as m^c € [1.6,2.0] 
GeV/c 2 and AE € [-0.14,0.14] GeV, excluding the 
blinded 3a ellipse. The toec vs. AE distributions 
of events inside the FB are modeled by 2-dimensional 
probability density functions (PDFs) summed over all 
background event types. The PDFs have correlations 
built in using Gaussian weights with an adaptive ker- 
nel estimation procedure [l8[. The shape of the Bhabha 
component is obtained using the data samples having 
cos# rcco ii < —0.8 from events selected in the e-tag sam- 
ple for the —> e ± 7 search, while the shapes of p + p~ , 
t + t~ and qq PDFs are obtained from their respective 
MC samples. 

The fractions of events for each background type are 
obtained from separate maximum likelihood fits to 41 
and 105 events inside the FB, respectively, for the two 
searches. We find (70 ± 15)% and (90 ± 8)% of the back- 
ground events are r-pair events. By integrating the to- 
tal PDF summed over background types only, we expect 
(1.6 ± 0.3) and (3.6 ± 0.4) events inside the 2a signal el- 
lipse for the two searches, where the quoted statistical 
errors are due to the sizes of the fitted dataset. 
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TABLE I: Means and resolutions of m^c and AE distributions for the signal MC events, the numbers of observed (obs) and 
expected (exp) events inside the 2a signal ellipse, the signal efficiencies (e), and the 90% C.L. upper limits (UL). 
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FIG. 1: The GSB and the 2a ellipse for t — ¥ e 7 (left) and r —¥ ^7 (right) decays in the mEC vs. AE plane. Data are 
shown as dots and contours containing 90% (50%) of signal MC events are shown as light- (dark-) shaded regions. 



As a cross-check, we integrate the total PDF over four 
2(j ellipses inside the FB, whose centers are shifted by 
±5(7 or by ±9tr along m,EC only. The numbers of ob- 
served events in each of these neighboring regions and 
their sums are consistent with the expected numbers of 
events, which are shown along with their statistical errors 
in Table HU 

TABLE II: Numbers of observed (obs) and expected (exp) 
numbers of background events along with statistical errors 
inside 2a ellipses whose centers are shifted by ±5<r and ±9cr 
in mEC only, and their sums. 
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14.9±0.8 



To obtain the systematic errors on the numbers of ex- 
pected background events, we fit the to.ec distributions 
of 32 and 81 data events inside the ±2<r band in AE over 
the GSB region but outside the blinded 3er ellipse. Vary- 
ing degrees of polynomial functions are used to model 



the Toec distributions, which are then integrated to ob- 
tain the number of expected events inside the 2a ellipse. 
The largest deviations between the predictions from 1- 
dimensional and 2-dimensional fits are used to set the to- 
tal uncertainties of 0.4 and 0.7 events on the background 
estimates. 

The systematic uncertainties in the signal selection and 
reconstruction efficiencies for — ¥ e ± 7 and r — ¥ /i 7 
decays due to the modeling of the variables entering the 
NN are 2.7% and 1.8%, respectively. Those due to the 
photon reconstruction efficiency are 1.8% for both de- 
cays, while those due to the signal-lepton track identi- 
fication are 2.3% and 2.7%, respectively. The contribu- 
tions due to the uncertainty in the signal-track momen- 
tum and signal-photon energy scale and resolution, es- 
timated by varying the peak position and resolution of 
the ttt-ec and AE distributions, are 6.4% and 6.2%, re- 
spectively. Other systematic uncertainties totaling less 
than 1.5% for both signal decay modes include those 
arising from trigger and filter efficiencies, tracking effi- 
ciencies, and the beam-energy scale and spread. We use 
approximately 10 6 MC events per channel, resulting in 
a negligible systematic uncertainty due to MC statistics. 
Although the signal MC has been modeled using a flat 
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phase space model, the efficiencies are insensitive to this 
assumption as demonstrated by considering the two ex- 
treme cases of V — A and V + A forms of interaction for 
the signal MC. All contributions to the systematic un- 
certainties are added in quadrature to give total relative 
systematic uncertainties on the efficiencies of 7.7% and 
7.4% for — > e ± 7 and t ± — > ^ ± 7 decays, respectively. 

We observe and 2 events for the ^ — > e ±r y and 
— > /i ± 7 searches inside the 2ct signal ellipse, respec- 
tively. As there is no evidence for a signal, we set a frc- 
quentist upper limit calculated using S^° L = N^° L /(N T e) 
to be 6(r± 6*7) < 3.3 x 1(T 8 and B(r ± ->■ ^7) 
< 4.4 x 10~ 8 at 90% C.L., where £ is the signal effi- 
ciency inside the 2a signal ellipse and Njj° L is the 90% 
C.L. upper limit on the number of signal events, esti- 
mated using the POLE program [191] . The upper limits 
which include all systematic uncertainties, are presented 
in Table HI along with signal efficiencies and numbers of 
observed and expected background events. These results 
supersede previous BaBAR results [3, US], reducing the 
upper limits by factors of 3.3 and 1.5, respectively, and 
are the most stringent limits on searches for lepton flavor 
violation in r — > e 7 and r — > /i 7 decays. 
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